The Fractional Quantum Hall States of Dirac Electrons in Graphene 
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We have investigated the fractional quantum Hall states of Dirac electrons in a graphene layer 
in different Landau levels. The relativistic nature of the energy dispersion relation of electrons in 
graphene significantly modifies the inter-electron interactions. This results in a specific dependence 
of the ground state energy and the energy gaps for electrons on the Landau-level index. For the 
valley-polarized states, i.e. ai v — 1/m, m being an odd integer, the energy gaps have the largest 
values in the n = 1 Landau level. For the valley-unpolarized states, e.g., for the 2/3 state, the 
energy gaps are suppressed for n = 1 as compared to those at n = 0. For both n — 1 and n — the 
ground state of the 2/3 system is fully valley-unpolarized. 
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A two-dimensional electron system in a single layer of 
graphite (graphene) is known to exhibit many remark- 
able properties. From the band structure studies 1] it 
was established early on that, to a good approximation, 
the energy dispersion of electrons in graphene is linear 
near the points at the corners of the Brillouin zone where 
the valence band and the conduction band meet. As a 
consequence, the low-energy excitations follow the Dirac- 
Weyl equations for massless relativistic particles 0. In 
an external magnetic field, the electron system also shows 
unique properties that are different from those of the 
standard non-relativistic electron systems |^ |^ Q . Re- 
cent experimental demonstration of some of those proper- 
ties, in particular, the discovery of the integer quantum 
Hall effect 1^ that was predicted in earlier theoreti- 
cal works H has caused intense interest in the electronic 
properties of the Dirac electrons in graphene j^j. How- 
ever, effects of electron correlations in this system have 
not been reported yet. In this paper, we report on the 
nature of the fractional quantum Hall states of Dirac elec- 
trons in graphene. 

A unit cell of the two-dimensional (2D) graphene hon- 
eycomb lattice contains two carbon atoms, say A and B. 
The dynamics of electrons in graphene is described by 
a tight-binding Hamiltonian with the nearest-neighbor 
hopping. In the continuum limit this Hamiltonian gen- 
erates the band structure with two tt bands and the 
Fermi levels are located at two inequivalent points, K = 
(27r/a)(l/3, 1/V3) and K' = (27r/a)(2/3, 0), of the first 
Brillouin zone, where a = 0.246 nm is the lattice con- 
stant. Near the points K and K' the electrons have a 
linear Dirac- Weyl ( "relativistic" ) type dispersion relation 
[112. Finally, in the continuum limit the electron wave 
function is described by the 8-component spinor, '^s,k,a, 
where s = ±1/2 is the spin index, k = K,K' is the valley 
index, and a — A, B is the sublattice index. Without 
the spin-orbit interaction fill IT^ the Hamiltonian is 
described by two 4x4 matrices for each component of 
the electron spin. In the presence of a magnetic field per- 



pendicular to the graphene plane the Hamiltonian matrix 
has the form 
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(1) 

where tt = p -I- eA/c, p is the two-dimensional momen- 
tum, A is the vector potential, and 7 is the band param- 
eter. The ordering used for the basis states in the non- 
interacting Hamiltonian is {K, A; K, B; K' , A; K' , B). 
The eigenf unctions of the Hamiltonian are specified by 
the Landau-level index n = 0, ±1, ±2, . . . and the intra- 
Landau level index m that is gauge dependent. Each 
Landau level is four-fold degenerate due to the spin and 
valley degrees of freedom. The corresponding wave func- 
tions for an electron in the two valleys K and K' are 
described by 



/ sgn(n)il"l V|n|-i \ 
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where C„ = 1 for n = and C„ = 1/V2 for n ^ 0. Here 
(j)n is the standard Landau wave function for a particle 
with non-relativistic parabolic dispersion relation in the 
n-th Landau level. From Eqs. 0-® it is clear that a 
specific feature of the relativistic dispersion law is the 
mixing of the non-relativistic Landau levels. This mix- 
ture is present only for n 7^ and strongly modifies the 
inter-electron interaction within a single Landau level. 
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In what follows, we study the partially occupied Lan- 
dau levels with fractional filling factors. Partial occu- 
pation of the Landau levels is realized by doping of the 
graphene layer. Experimentally, different filling factors 
of the Landau levels are achieved by varying the applied 
magnetic field at a fixed electron concentration. In this 
case the ground state of the system and the excitation 
spectrum are fully determined by the inter-electron inter- 
actions. For the non-relativistic case this results in the 
incompressible fractional quantum Hall effect (FQHE) 
states at the fractional filling factors [11 m. Proper- 
ties of these states are completely described by Haldane's 
pseudopotentials, Vm which are the energies of two 
electrons with relative angular momentum m. The pseu- 
dopotentials for the n-th Landau level can be presented 
as 



^qV{q)[F^{q)fL„.{q^)e-^\ (4) 



where Lm{x) are the Laguerre polynomials, V{q) = 
2^6^ / (nlq) is the Coulomb interaction in the momen- 
tum space, K is the dielectric constant, / is the magnetic 
length, and Fn{q) is the form factor corresponding to 
the n-th Landau level. For relativistic electrons the form 
factor is given by the expression 
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while for the non-relativistic particles the form factors 
in Eq. ^ are Fn{q) = Ln (9^/2). This means that the 
inter-electron interactions for the relativistic and non- 
relativistic electrons are the same for n = and are dif- 
ferent for n > 0. In what follows, all energies are 
expressed in units of the Coulomb energy Sc = / kI. 

In Fig. 1 we compare the pseudopotentials calculated 
from Eq. Q for the relativistic and the non-relativistic 
cases. For the relativistic electrons [Fig. 1 (a)] we notice 
a clear suppression of the pseudopotential for n = 1 as 
compared to that at n = 0, only for m — Q, i.e. when 
both electrons are at the same spatial point. For all other 
m, we have the inequality Vm'' > Vm\ This is differ- 
ent from the non-relativistic case where the pseudopo- 
tential is suppressed also for m — 1, i.e. v}^^ < Vi^\ 
We also see in Fig. 1 (a) that although the relativistic 
wave functions at n = 1 is the "mixture" of the n = 
and n — 1 non-relativistic wave functions, the relativistic 
pseudopotential is not the average of the corresponding 
non-relativistic pseudopotentials. This is clearly seen for 
TO = 1 where the relativistic pseudopotential at n = 1 is 
larger than the non-relativistic one for both n = and 
n = 1. In Fig. 1 (b) the relativistic pseudopotentials 
are shown for different Landau levels. Here the special 
case is TO = 1 where the dependence of the pseudopo- 
tential on the Landau-level index is non-monotonic, viz., 
the pseudopotential has the maximum value ai n = 1. 



,n=0 (relativistic) 

. n=1 (relativistic) 
. n=1 (non-relativistic) 
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FIG. 1: The pseudopotentials [Eq. (0J] as a function of the 
relative angular momentum (a) for relativistic and for non- 
relativistic 2D electrons in the first two Landau levels, and 
(b) for relativistic electrons in various Landau levels. 



At all the other m values the trend is the same as for 
the non-relativistic case, i.e., for m — the pseudopo- 
tential decreases with increasing n, while at m > 1 the 
pseudopotentials increase with n. 

With the pseudopotentials for Dirac electrons at hand, 
we now evaluate the energy spectra of the many-electron 
states at fractional fillings of the Landau level. The cal- 
culations have been done in the spherical geometry [l^ 
with the pseudopotentials given by Eq. Q). The radius of 
the sphere R is related to 25* of magnetic fluxes through 
the sphere in units of the flux quanta as i? = y/Sl. The 
single-electron states are characterized by the angular 
momentum S, and its z component Sz- For a given num- 
ber of electrons N, the parameter S determines the filling 
factor of the Landau level. Due to the spherical symme- 
try of the problem, the many-particle states are described 
by the total angular momentum L, and its z component, 
while the energy depends only on L. At first we study the 
system with the fractional filling factor i' = 1/3, which 
corresponds to the 1/3 FQHE. In the spherical geometry, 
the 1/3-FQHE state is realized at 5* = (3/2)(Af - 1). If 
the electron system is fully spin and valley polarized then 
we should expect the ground state to be in the Laugh- 
lin state jl^ which is separated from the excited states 
by a finite gap. We calculated the energy spectra of a 
finite-size system by finding the lowest eigenvalues and 
eigenvectors of the interaction Hamiltonian matrix ^1§\. 
In these calculations we take into account the interaction 
between all the electrons of the partially occupied Lan- 
dau levels. We have also addressed the question of polar- 
ization of the many-particle state. We assume that at a 
high magnetic field the system is always spin-polarized. 
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FIG. 3: Energy spectra of the eight-electron u = 2/3-FQHE 
system, shown for different Landau levels: n — (stars) and 
n = 1 (dots). The system is valley-unpolarized and fully spin- 
polarized. The flux quanta is 2S = 11. 



FIG. 2: Energy spectra of the eight-electron u — l/S-FQHE 
system, shown for different Landau levels: (a.) n — (stars) 
and n — 1 (dots), and (b) n — 2. The system is fully spin and 
valley polarized. The flux quanta is 2S = 21. 

In this case the system can be either valley-polarized or 
valley-unpolarized. Similar to the spin polarization of 
the standard FQHE states the valley polarization of the 
graphene system depends on the inter-valley asymmetry 
(same as the Zeeman splitting for the non-relativistic sys- 
tem) and filling factor of the Landau level. Namely, if the 
inter-valley splitting due to the inter-valley asymmetry 
is large then the system is always valley-polarized, while 
otherwise the valley polarization of the system depends 
on the filling factor of the Landau level. 

In Fig. 2 we show the calculated energy spectra for the 
l/S-FQHE state at different Landau levels. Here, the 
1/3-FQHE state at the n-th Landau level is defined as 
the state corresponding to the 1/3 filling factor (single 
valley, single spin) of the n-th Landau level, while all 
the lower energy Landau levels are completely occupied. 

Since the relativistic pseudopotential Vm "' for n = is 
the same as for the non-relativistic one, the 1/3 state 
and the corresponding energy gap (in units of £c) will 
be the same in both cases. The deviation from the non- 
relativistic system occurs only at higher Landau levels. 
In Fig. 2 (a) the energy gap of the 1/3-state at n = 1 is 
noticeably enhanced compared to that at n = 0. This is 
a direct manifestation of the specific dependence of pseu- 

(n) 

dopotenials Vm on the Landau-level index. Due to the 
asymmetry of the electron wave functions, the spectra of 
the 1/3-FQHE state is mainly determined by the relative 
value of v/"''' and V^"'^ pseudopotentials, which have the 
highest value at n = 1. The energy spectra of the 1/3- 
FQHE state at n — 2, shown in Fig. 2 (b), demonstrate 
a strong suppession of the gap when compared to the 
n = 1 and n = FQHE states. We therefore conclude 



that the 1/3-FQHE state in graphene is most stable at 
n = I. Hence the inter-electron interaction effects are 
more pronounced at n = 1. Interestingly, this tendency 
is just the opposite to that of the non-relativistic system, 
where the excitation gap decreases monotonically with 
incresing Landau-level index {14] . 

The results in Fig. 2 correspond to a completely 
spin and valley polarized system. This polarization is 
achieved at a high magnetic field due to the Zeeman split- 
ting and the valley asymmetry. The inter-valley asym- 
metry has two sources: the first one is due to interaction- 
induced "backscattering" between different valleys 
while the second one is due to the asymmetry in the 
lattice-scale interactions within the two sublattices of 
graphene Since the positions of the electrons in 

two sublattices are shifted, the interaction between the 
electrons in the different sublattices is weaker than the 
interaction between the electrons in the same sublattice. 
Both effects vary as {a /I), so they become more relevant 
at higher magnetic fields or at a smaller magnetic length. 

The same picture holds for the other FQHE states of 
the type 1/m, i.e., the state is most stable at n = 1. The 
new aspects of the interaction physics arise at other filling 
factors as well, when the lowest energy states are spin- 
unpolarized for the non-relativistic single- valley systems. 
The simplest example is v = 2/3. In this case the ground 
state of the non-relativistic electrons is spin-unpolarized 
at a small Zeeman splitting. The transition from the 
spin-polarized to the spin-unpolarized ground states in a 
tilted magnetic field is well established both experimen- 
tally and theoretically for non-relativistic electrons 15]. 

Here we address the problem of the polarized and un- 
polarized states in the graphene system. We assume 
that the Zeeman splitting is large enough so that all 
the states are spin-polarized. At the same time the 
inter-valley asymmetry is small and electrons can occupy 
both the valleys. In this case the valley-polarized and 
valley-upolarized 2/3-FQHE states become relevant. For 
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a valley-polarized system the excitation spectra and the 
ground state properties of the 2/3-state are the same as 
those for v — 1/i due to the particle- hole symmetry. Sim- 
ilar to the 1/3 case we obtain numerically an enhance- 
ment of the excitation gap at n = 1 as compared to that 
at n = 0. The more complicated situation occurs for the 
valley-unpolarized system. First we compare, just as for 
the non-relativistic system T^l , the energy of the ground 
states of the polarized and the unpolarized systems. In 
the spherical geometry the polarized and the unpolarized 
states are realized for different sizes of the sphere, i.e., for 
different flux quanta through the sphere. For the 2/3- 
FQHE system the polarized state occurs at 2S = 3N/2, 
while the unpolarized one is at 2S — 3N/2 — 1. Due 
to different size of the sphere in these two systems the 
finite-size corrections to the magnetic length should be 
introduced, I' = {v2S/NY/H ^Q^. 

We have calculated the ground state energies for the 
valley-polarized and the valley-unpolarized graphene sys- 
tem in the n = 1 and the n = Landau levels for a eight- 
electron system in the spherical geometry. We found that 
for both n = and n = 1 the valley-unpolarized state 
has the lower energy. The ground state energy per par- 
ticle in the unpolarized system is lower than that for the 
polarized system by 0.073ec for n = and by 0.053£c for 
n = 1. Therefore, the valley-unpolarized state is more 
favorable for n = 0. Here the gap between the polar- 
ized and the unpolarized states is suppressed for n = 1 
as compared to the n = case. This is opposite to the 
completely polarized system (Fig. 2), where the effects 
of interaction is the strongest for n = 1. Suppression of 
the interaction effects in an unpolarized system for n — 1 
is also illustrated in Fig. 3, where the excitation spec- 
tra of the valley-unpolarized system is shown for n = 
and n = 1. A strong suppression of the excitation gaps 
for n = 1 is clearly visible here. The origin of this sup- 



pression can be understood from the dependence of the 
pseudopotentials kI"' on the relative angular momentum 
TO in different Landau levels. Due to the Pauli exclusion 
principle the energetic properties of the polarized state 
is determined only by the pseudopotentials with odd an- 
gular momenta, m = 1,3,5,.... These pseudopotentials 
have the largest values for n = 1, which results in the 
strongest interaction effects for n = 1, in a polarized sys- 
tem. For an unpolarized system the properties of the 
ground and excited states depend on all the pseudopo- 

(n) 

tentials. Since Vm at to = is strongly suppressed for 
n = 1 as compared to the n = case, we expect a sup- 
pression of interaction effects in an unpolarized system 
in the n = I Landau level. 

In conclusion, the relativistic nature of the energy dis- 
persion of electrons in the graphene plane modifies the 
inter-electron interactions significantly. This results in a 
unique dependence of the ground state energy and the 
energy gaps of the graphene systems on the Landau-level 
index. For the valley-polarized states, i.e. a.i v = 1/m, 
the FQHE gaps have the largest values for n = 1. Based 
on these studies we conclude that the FQHE at = 1/to 
should be observed experimentally for both n = and 
n = 1, perhaps in a higher mobility system. For the 
valley-unpolarized states, e.g., for the 2/3-FQHE state, 
the energy gaps are suppressed at n = 1 as compared 
to that for the n = level. For both n = 1 and 
n — the ground state of the 2/3-FQHE system is fully 
unpolarized. The inter-valley asymmetry will result in 
transitions between the valley-polarized and the valley- 
unpolarized states. 
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